Europhysics Letters PREPRINT 



■^ ' Multi-terminal transport through a quantum dot in the 
Q ■ Coulomb blockade regime 

^ ■ R. Leturcq^(*), D. Graf^, T. Ihn^, K. Ensslin\ D. D. Driscoll^ and A. C. 

^ I GOSSARD^ 

^ Solid State Physics Laboratory - ETH Zurich, 8093 Ziirich, Switzerland 
CN , 2 Materials Department, University of California - Santa Barbara, Ca 93106, USA 



C^ ' PACS. 73.23.Hk - Coulomb blockade; single-electron tunneling. 

^ I PACS. 73.63.Kv - Quantum dots. 

Abstract. - Three terminal tunnelling experiments on quantum dots in the Coulomb-blockade 

j^ , regime allow a quantitative determination of the coupling strength of individual quantum states 

to the leads. Exploiting this insight we have observed independent fluctuations of the coupling 

strengths as a function of electron number and magnetic field due to changes in the shape of the 

'■^ , wave function in the dot. Such a detailed understanding and control of the dot-lead coupling 

$Ih ' can be extended to more complex systems such as coupled dots, and is essential for building 

O I functional quantum electronic systems. 



^ ' In a standard two-terminal experiment with a single quantum dot in the Coulomb blockade 

^V , regime [1], the current in a conductance resonance is determined by the average coupling of 
^— ^ ' the electron wave function in the dot with the corresponding wave functions in both leads. In 
\^ , the linear regime, such an experiment does not allow to determine the individual coupling of 
^D ' the wave function in the dot to each terminal. Here we demonstrate that in the single-level 
^r , tunnelling regime of the Coulomb blockade it is possible to deduce the individual coupling 
strengths from the dot to the leads if three or more terminals are connected to the dot. It 
ti . is possible to determine the conductance matrix of the quantum dot, and to calculate the 
individual tunnelling rates from the dot to each lead. For weak coupling, the magnitude 
of the tunnelling rates of a given terminal is found to vary independently of the two other 
'^ ' tunnelling rates when the number of electrons in the dot is changed. This result can be related 
G , to the chaotic nature of the wave function in the dot. The fluctuations of the shape of the 
wave function in the dot due to quantum interference is directly observed via the magnetic 
field dependence of the coupling strengths. Finally, level broadening beyond thermal effects 
becomes visible as the coupling strength of a single lead increases. This set-up allows then to 
tune the tunnelling rates from the dot into the three terminals individually. 
^ ' The sample has been fabricated on an AlGaAs-GaAs heterostructure containing a two- 

dimensional electron gas (2DEG) 34 nm below the sample surface. A back gate situated 1.4 
/xm below the 2DEG allows to tune the electron density. All measurements presented here 
were performed at a back gate voltage of -1.4 V, giving a 2DEG density of 3.7 x 10^^ cm^^ and 
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a mobility of 200'000 cm^/Vs at T — 4.2 K. The surface of the heterostructure was locaUy 
oxidised by applying a voltage between the conductive tip of an atomic force microscope 
(AFM) and the 2DEG [2]. The electron gas is depleted below the oxide lines. This patterning 
technique was used in other studies for defining high-quality quantum dots [3,4]. The details 
of the fabrication process, which are crucial for the high electronic quality of the quantum dot, 
are described in ref. [5]. Figurena) shows the oxide lines defining the quantum dot. The AFM 
image was taken with an unbiased tip directly after the oxidation process. The width of the 
four quantum point contacts connecting the dot to the four reservoirs numbered 1 through 4 
is controlled by voltages applied to the lateral gate electrodes LGl through LG4. The number 
of electrons in the dot can be tuned via the lateral plunger gate PG. In this paper we focus on 
the dot being connected to three terminals. The point contact connecting the dot to reservoir 
4 is completely closed, and gate LG4 can be used for controlling the number of electrons in 
the dot. Measurements of Coulomb diamonds reveal a charging energy Ec ~ 0.5 meV and 
an average single-particle level spacing A w 35 /ieV, compatible with an electronic dot area 
of 400 nm x 250 nm estimated from the lithography pattern and the lateral depletion. An 
electronic temperature of 90 mK is deduced from the width of the Coulomb peaks, as it will be 
explained later. The quantum dot can be tuned into the quantum Coulomb blockade regime 
with the mean single-particle level spacing A being much larger than the thermal energy ksT 
and the level broadening hT. 

Figure ^h) shows the measurement set-up. A dc bias voltage of 10 /iV is applied to 
one terminal of the dot {e.g. Vbiasi), while the two other terminals are grounded {e.g. 
Vbias2 = Vbiass — 0). Current- voltagc converters are used to measure the currents through 
each terminal. We present the experimental results in the following way: for a bias voltage 
applied to one terminal, the currents through all three terminals are measured. In order to 
minimise the influence of possible offsets, measurements for positive and negative bias are 
averaged. Then, by applying the bias successively to each of the three terminals, we obtain 
nine different current measurements. In linear response theory, these nine currents correspond 
to the nine elements of the conductance matrix G of the three-terminal system: 

Gil Gi2 Gi3 

= ( G21 G22 G23 

G31 G32 G33 

In a preliminary experiment on a strongly coupled dot some of the Gy 's have been measured 

[6]. 

Figuresnc)-e) show the nine conductances Gy as a function of the gate voltage Vlg4 that 
controls the number of electrons on the dot (in each panel, one curve is laterally offset by 
3 mV for clarity). The positions of all corresponding Coulomb resonances agree within less 
than 1/10 of the peak width {i.e. less than 5 fJ,eV), indicating that the same energy level in 
the dot is probed in all configurations. Current conservation implies X]i=i ^y ~ ^ ^^r all j. 
Additionally, if all voltages are set to the same value no current should fiow: J2i=i ^ij = 
for all i. Figures ^f) a-nd g) show that these two sum rules are obeyed by the experimental 
data with a relative accuracy better than 10% of the highest current level. 

In the case of very low temperature and weak coupling, one can use the theory for lowest 
order sequential tunnelling including interaction effects in the dot, generalised to the case of 
more than two terminals. Following Beenakker [7], we find from the rate equation approach 
that, in the regime of weak coupling {hT <^ ksT) and in the single-level transport regime 
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Fig. 1 - a) Micrograph of the four-terminal quantum dot, the black lines being the oxide lines. The 
four leads (labelled 1 to 4) to the dot can be tuned through the four lateral gates LGl to LG4. The 
plunger gate PG tunes the number of electrons in the dot. b) Measurement set-up using three of 
the four terminals. The quantum dot (QD) is connected to the three leads through tunnel barriers 
with tunnelling rates Fi to F3. c)-e) Measurement of Coulomb blockade resonances in the three 
configurations, when applying a bias voltage to one lead, the others being grounded. For each plot, 
the negative conductance corresponds to the current in the biased lead (1 in c), 2 in d) and 3 in e)), 
and the positive conductances to the two currents in the grounded leads. One of them is laterally 
offset by a constant of +3 mV to make the presentation more transparent, f) Current sum rule for 
the conductance at the resonances presented in panels c) to e), obtained by adding the terms of each 
column of the conductance matrix: Gn + G21 + G31 (+), G12 + G22 + G32 (o) and G13 + G23 + G33 
(A), g) Voltage sum rule obtained by adding the terms of each row of the conductance matrix: 
Gil + G12 -f Gi3 (-H), G21 + G22 + G23 (o) and G31 + G32 + G33 (A). 



(fc^T ^ A), the elements of the conductance matrix G are given by 
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with Ffe being the tunneUing rate from the dot to lead k (see fig-Hb)), and S — eaLG4( Vlg4, res — 
^LG4)i with ai(34 the lever arm of gate LG4, determined by the measurement of Coulomb 
diamonds. Each peak shown in fig. n]c)-e) is fitted with eq. J^l or |(2Jl in order to deduce its 
position, its maximum and its width. From the maxima of the peaks, we calculate the indi- 
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Fig. 2 - Tunnelling rates hT^ for individual resonances for lead 1 (+), lead 2 (o) and lead 3 (A), and 
corresponding conductances Gk (right hand axes). Inset: Peak width at half maximum (FWHM) for 
the resonance measured in Gn (+), G22 (o), G33 (A), and mean peak width of the nine conductance 
matrix elements (thick line), a) For low Vlga, i-e. weak coupling, b) At higher Vlga, i-e. intermediate 
coupling. Inset: The dashed thick line at 27 /ieV corresponds to a thermally broadened peak at 90 
mK, the temperature corresponding to the mean peak width for a more closed dot (inset of panel a)). 



vidual tunnelling rates Ffc. The tunnelling rates depend on the overlap of the wave function 
in the dot with the wave function in lead k. This overlap depends on the shape of the wave 
function in the dot, and therefore on the interference pattern in the dot. 

The tunnelling rates measured at resonance in the weak coupling regime as a function of 
the gate voltage Vlg4 are shown in fig. [21a). Such an experimental determination of tunnelling 
rates is unique to a three (or more) terminal quantum dot and not possible in conventional 
two-terminal experiments in the linear regime. The values of the tunnelling rates fluctuate 
strongly and independently [8] , also when the mean coupling strengths are similar for all leads 
(fig. [21a)). Only for more positive gate voltages (fig. [21b)), where the dot starts to open because 
of electrostatic coupling of the plunger gate to the point contacts, the tunnelling rates start 
to increase systematically. However, in this case, the single-level tunnelling regime might not 
be valid anymore because of the higher tunnelling rates, and the determination of individual 
tunnelling rates is no longer valid due to co-tunnelling. 

It is known that states in a dot change also as a function of magnetic field due to changes 
in the quantum interference pattern. We have measured the conductance matrix around one 
resonance as a function of the magnetic field B and the gate voltage Vlg4, as shown in fig. [31 
The fact that the peak positions fluctuate the same way as a function of the magnetic fleld 
for all configurations confirms that the three leads couple to the same state. Following each 
peak, we have extracted the peak heights and plot the tunnelling rates as a function of the 
magnetic field in fig. 21 Again, the three tunnelling rates fiuctuate independently. 

We like to note that, at a resonance, the G^ defined in eqs. (^ and ^ are formally 
equivalent to the conductance matrix elements of a classical star-shaped conductance network, 
the three conductances connecting each lead k to the center of the dot being defined as 
Gk = {e^ /h){hTk)/{^kBT). We can therefore represent hVu in figs. [21 and ^ as equivalent 
conductances Gk (right hand axes). 

Does coherent tunnelling change this analysis? It has been shown experimentally that 
tunnelling through a quantum dot is at least partially coherent [9]. Coherence may have an 
influence on the dot conductance if interference is possible [10]. If a particular level spacing 
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Fig. 3 - Absolute value of the conductance matrix elements dj vs. Vlg^ and magnetic field B around 
a resonance. The conductance is represented with a linear gray scale with white corresponding to 
and black to > 0.025 e^ /h. 

Fig. 4 - Tunnelling rates at a resonance hVk, for lead 1 (+), lead 2 (o) and lead 3 (A), vs. magnetic 
field B, and corresponding conductances Gk (right hand axes). The magnetic field necessary to add 
a flux quantum through the dot area is 40 mT. Inset: Peak width at half maximum (FWHM) for the 
resonance measured on Gn (+), G22 (o), G33 (A), and mean of the nine conductances (thick line). 



5e becomes significantly smaller than the average level spacing A « 35 /ieV, the condition 
hT > 5e leads to level mixing due to interference. For such a case our relation between 
the Gij and the Ffe has to be modified. However, if the coupling is weak, hT ^ A, no 
interference is possible. In this case, it has been shown at least for the two-terminal experiment 
that the expression for coherent tunnelling reduces to the known equations for sequential 
tunnelHng [11]. 

Equation. (Q implies that Gij = Gji because Gij oc FiFj. This agrees with the gener- 
alized Onsager relations Gij{B) — Gji{—B) at zero magnetic field. At finite magnetic field 
this symmetry can not be observed in general [12], an extreme example being an electron 
focusing experiment [13], and eq. |^ may not have the most general form required for finite 
B. However, we find that in our experiment Gij{B) — Gji{B) for magnetic fields below 0.3 T 
(see fig. 13, which implies Gij{B) ~ Gij{—B) through the Onsager relations. In the absence 
of a more appropriate theory we therefore employ eq. (Q empirically also at finite B. Never- 
theless we may ask whether the symmetry of the conductance matrix in our measurement is 
a general property of any three-terminal quantum dot or the result of particular microscopic 
properties of our system. We speculate that the observed symmetry is a general property of 
a multi-terminal system coupled weakly to the leads, because the tunnelling matrix elements 
depending on the overlap of dot and lead wave functions are even in magnetic field. In lowest 
order in the tunnelling the Gij will therefore be symmetric in B. However, we believe that 
this question has to be settled eventually with new experiments and a thorough theoretical 
analysis. 

From the measurement of individual dot-lead coupling strengths, we can address three 
questions. How can we probe changes in the wave function in the dot? How is the level 
broadening affected by opening a single contact? And how can the coupling strengths be 
monitored on a quantitative level? 

Peak height fluctuations attributed to fluctuations of the shape of quasi-bound states in 
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chaotic dots [1] have been extensively studied in two-terminal quantum dots [14,15]. Cal- 
culations based on random matrix theory are in reasonable agreement with experimental re- 
sults [10, 14, 16]. However, two-terminal experiments can only give information on the global 
conductance of the entire system. In a three-terminal setup this picture can be probed more 
directly by looking at the spatial distribution of the wave function. In our experiment, these 
spatial fluctuations are observed as uncorrelated tunnelling rates because the distance between 
leads (« 400 nm) is much larger than the Fermi wave length (w 40 nm) [11]. A magnetic field 
perpendicular to the 2DEG changes the shape of the wave function, as it has been suggested 
for the interpretation of the peak height fluctuations as a function of the magnetic held in 
two-terminal experiments [14]. In our experiment, independent fluctuations of the individual 
tunnelling rates on a scale of 40 mT, corresponding to one flux quantum added within the dot 
area, are a direct consequence of the spatial fluctuations of the wave function due to changes 
of the interference pattern in the dot. 

In order to check the influence of the coupling to the leads on the level broadening, an effect 
that is beyond eqs. ^ and Q, we have also carried out the analysis of the peak width for all 
configurations. For weak coupling, hT <§; fc^T -C A <C Eq-, the peak width is approximately 
constant [inset of fig.|21a)]. The full width at half maximum (FWHM) of 27 /LteV corresponds 
to the width of temperature broadened peaks at 90 mK [1]. At larger gate voltage [inset 
of fig. Elb)] , the coupling to the leads becomes stronger and the dot enters an intermediate 
coupling regime {hT < ksT). The peak width increases continuously as a function of gate 
voltage due to an increase of the level broadening. This shows directly the influence of the 
coupling strength on the peak width. Although the tunnelling rate of lead 1 increases more 
than the others [flg. Elb)] , the widths of all Coulomb peaks increase the same way. This means 
that the level broadening due to the coupling to one lead can be seen in the width of all the 
other conductances, as expected if all leads are coupled to the same state. Here we point 
out that this last result only holds on a qualitative level, since the dot is no longer in the 
single-level tunnelling regime. In the inset of fig. Old), the same analysis is carried out for the 
magnetic field dependence. The peak width does not change when increasing the magnetic 
field, although the tunnelling rates vary strongly, meaning that the dot is still in the weak 
coupling regime. 

The implementation of spin-based quantum information processing [17-19] with quantum 
dots requires more complex and coupled quantum systems. Two coupled quantum dots em- 
bedded in an Aharonov-Bohm interferometer are suggested for the detection of entanglement 
of spins states [18]. A three-terminal dot could be a step towards building an electronic en- 
tanglcr [20]. For such devices, the quantitative control of individual coupling strengths will 
be a prerequisite for the desired functionality. We have demonstrated that a dot in the weak 
coupling regime (see fig. Ela)), where tunnelling through a single level is expected to occur, 
can be tuned into a regime where the average coupling strengths to all leads are approximately 
the same. However, fiuctuations of the shape of the wave function in the dot do not allow 
to control each coupling strength independently. Quantum rings, which have a more regular 
energy level spectrum rather than a chaotic one, even for a large number of electrons [4], could 
be better candidates. In general, our method of measuring multi-terminal quantum dots pro- 
vides a means to deduce the coupling strengths in more complex systems. It is also suitable for 
studying interference effects on the peak height statistics in chaotic quantum dots [10, 14, 15] 
by working directly with the individual tunnelling rates rather than with the two-terminal 
conductance. 
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